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SUMI-IARY 


An e>j)erimental investigation of the general instability of reinforced 
thin-walled metal cylinders was carried out at the California Institute of 
Technology. The basic parameters involved were the spacing and sectional 
properties of the stiffening elements^ the wall thickness, and the diame- 
ter of the cylinder. An analysis of the experimental data led to a suita- 
ble parameter for estimating the general— instability stress of reinforced 
metal cylinders when subjected to pure torsion loading. 


INTRODUCTION 


The present report deals with the experimental investigation of the 
general instability of metal cylinders subjected to xnire torsion loading. 
Reports on other loading conditions, that is, pure bending, combined bend- 
ing and transverse shear, and combined bending and torsion, have been pub- 
lished previously (references 1 to 7)* . 

Inasmuch as a condition of a, pure torsional load seldom arises in the 
design of fuselage or wing structures, the problem of pure torsion as such 
might not warrant ah investigation. However, under a combined loading of 
bending and torsion (reference 7) the ultimate load of the stiffened cyl- 
inder is dependent on the ratio of the shearing stress at failure for com- 
bined loading to the shearing stress at failure for pure torsion. Hence, 
in order to predict the ultimate strength of a stiffened metal cylinder 
subjected to combined bending and torsion, a knowledge of the ultimate 
strength of the cylinder when subjected to. a pure torsion loading is nec— • 
essary. 

Because of the nonlinearity of the buckling problem of stiffened cyl- 
inders (cf. references 5 ®nd 8), no attempt has been made to give a theo- 
retical treatment of the problem. As given in the present repojrb, the 
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paraiaeter for predicting the ultimate strength of stiffened metal cylin 
ders subjected to torsion loads is based on ah analysis of the experi- 
mental results and on the existing theory of unstiffened metal cylinders. 
This method was preferred over that of a linearized theory which cannot 
correctly describe the behavior of the structure. The results of a lin- 
ear theory would have to be modified and corrected to bring it Into 
agreement with the experimental observations and thus the theory would 
immediately be rendered an empirical method. 

This investigation was carried out by the California Institute of 
Technology under the sponsorship and with the financial assistance of 
the National Advisory Committee for Aeronautics. 


SYMBOLS 

Mj applied torsional moment, inch— pounds 

applied torsional moment at skin buckling, inch— pounds 

^max torsional moment at failure, inch-pounds 

T shearii^g stress in the sheet covering, pounds per square inch 

T-_ buckling shear stress of the sheet covering, pounds per square inch 

cr 

''’max shearing stress in the sheet covering at failure, pounds per square 
inch ■ ' 

tensile stress of diagonal— tension field, potmds per square inch 

V • • ... 

P tension load resulting from ct acting over a unit circximferential 
width of sheet; acts in the direction of cos pounds 

• ^ ! ‘ . . ..i* 

Pg longitudinal con 5 )onent of P, pounds - • ; 

Py circumferential component of P, pounds • • 

a angle of diagonal— tension field with reference to a longitudinal 
stiffener 

strain in a longitudinal stiffener 

G shear modulus, pounds per square inch 
E Young’s modulus, pounds per square inch 
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Gg effective shear modulus, pounds per square inch 

t thickness of sheet covering, inches 

A area enclosed by sheet covering, square inches 
E radius of cylinder, inches 

L length of cylinder, inches . , ; 

b . spacing of longitudinals, inches 
d spacing of frames, inches 

radius of gyration of a longitudinal and effective sheet, inches 
p radius of gyration of a frame and effective sheet, inches 


DESCEIPTION OF APPARATUS 

All tests were conducted in the combined bending and torsion machine 
of the structures laboratory' at GALCIT. This machine has a maximum 
capacity of 500,000 inch— pounds. A detailed description of the machine 
is given in reference 2. A photograph (fig. 1) shows a 10-inch-diameter 
cylinder mounted for a pure— torsion test. 

The wire-strain-gage equipment, as used for the strain measurements, 
has been described in reference 6. The circumferential position of the 
various gages is shown in figure 2. Longitudinally the gage is mounted 
at the center of the specimen. 


TEST PROCEDURE 

The test specimens were all circular reinforced metal cylinders. 

Two sizes of specimens were tested, one series having a diameter of 32 
inches and a length of 64 inches, and the other series, a diameter of 20 
inches and a length of 40 inches. The sheet covering was 0.010-, 0.015—, 
and 0.020— inch-thick 24ST dural. The longitudinal reinforcing members 
consisted of round 24ST dural tubing drawn to an elliptical shape and 
the frames consisted of rectangular bars of 24ST Alclad. Longitudinals 
of three different wall thicknesses and frames of two different sizes 
were used. A sketch of the reinforcing members is presented in figure 3* 
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A diagram showing a shoet panel bounded by two longitudinal stringers and 
end mounting rings is presented in figure 4(a). The relationship of the 
forces acting on an element of sheet of unit width is given in figure 4(b). 
In table I the sheet thicloiess, the reinforcing members, and their spaciug 
are listed for each specimen. Stress— strain curves of the materials used 
in the tests are shown in figures 5 ho 12. 

All test specimens were tested in the combined torsion and bending 
machine mentioned previously . In the pure torsion loading, the bending 
arms and one torsion arm were locked in place, the load being applied to 
the second torsion arm. In all tests, shear stresses in the sheet cover- 
ing were not measured, inasmuch as it was assumed that a calculation of 
the shear stresses based on the classical formula 

T = Mr/2At (1) 

would be sufficiently acciirate. Strain measurements of the induced 
stresses in the longitudinal members were made on all specimens, the 
strain measurements being made by means of electric strain gages as de- 
scribed in reference 6. For a number of specimens, the unit angular 
deflection as a function of the applied torque was measured. By measu3>- 
ing the differential displacement between two pointers mounted on the 
cylinder, the unit angular deflection could be calculated. The pointers 
consisted of triangular frsmie structures and are shown mounted on the 
specimen in figure 13- Mounting these pointers on the cylinder just in- 
board of the end-rings ensured that the entire measured deformation oc- 
curred in the cylinder. 

The effective shear modulus was calculated in the folijowing manner: 

If 

I’ length of the pointer measured from the center of the cylinder, 
inches 

6 measured differential displacement between the two pointer, inches 
I distance between the pointers, inches 
then the unit angular displacement is given by 

e ^ 

IV 


and the effective shear modulus is 
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( 2 ) 


Some attempts were made to measure the wave fonn of the buckled 
sheet as a function of the applied torque; however, considerable diffi- 
culty was encountered in these measurements , primarily because of the 
nonuniform buckling which occurred over the cylinder. It was found that, 
in the testing of a cylinder, buckling of individual panels, that is, a 
sheet panel bounded by two frames and two longitudinals, would vary by as 
much as 100 percent . For this reason it was thought that a measurement 
of the wave form and wave amplitude as a function of the applied load 
would be rather meaningless. Photographs of a number of typical failures 
are shown in figures 13 to 17 . 


TEST RESULTS 


Results of the induced longitudinal strain as a function of the ap- 
plied torsional moment are shown for 17 specimens in figures I 8 to 
Each curve represents the average of two strain gages mounted on diamet- 
rically opposite sides of the longitudinal stiffener. 

The induced strain is associated with the shear buckles in the sheet 
covering. It can be readily seen that upon buckling of the sheet, longi- 
tudinal forces which act on the end mounting rings are introduced by the 
sheet covering. These forces have a tendency to pull the end mounting 
rings together and, since such a motion is resisted by the longitudinal 
stiffeners, the resvilt is an induced strain in these members. In the 
appendix some calculations have been worked out on the magnitude of the 
resulting strains . The calculated values do not agree too well numerir- 
cally with the measured values; however, the order of magnitude is cor- 
rect. 

The results of induced strain as a function of applied torque have 
been cross— plotted as shown in figures 35 tp 46. The strain was plotted 
radially outward with the cylinder circumferences as the zero reference 
line. These curves indicate that in general the induced strain is dis- 
tributed uniformly around the cylinder. 

In addition to the foregoing test data, the measurements of the ef- 
fective shear modulus as a function of the shearing stress in the sheet 
covering are given in figures 47 to 53* At the lower slear-stress values, 
the results are not very reliable because of the difficulty in measuring 
the resulting small deflections . At the higher values of t the results 
are sufficiently accurate to indicate the order of magnitude of the 
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effective shear modulus which may result in an actual structure. It is 
of interest to note that the effective shear modulus Gj; is considerahly 
less than 0.625G, which is the theoretical value for the fully developed 
diagonal— tension field of a flat plate. 

In presenting these experimental data, it would be more apporpriate 
to use the ratio rather than T or Mj since all the measured 

quantities are associated with buckling. However, as previously men- 
tioned, for the cylinders tested, buckling. did not occur uniformly over 
the cylinder at a definite load. The buckling process was progressive 
and uniform buckling occurred, for some specimens, at more than twice the 
load at which the first buckles appeared. For this reason, it was not 
possible to establish a definite observed buckling value. 


In detenoining a parameter for predicting the general instability 
failure in torsion, the same general procedure was followed as was used 
in the pure— bending Investigation of reference 5* The variables to be 
considered are the same as those of the pure— bending problem and can 
again be divided into two groups, namely, those dealing with the geoia- 
etry of the structure and those involving the sectional properties of 
the stiffening elements as well as the sheet covering. The geometrical 
variables are the longitudinal spacing b, the frame spacing d, the 
diameter, and the length of the cylinder. The second group of- variables 
includes the section properties of the longitudinals and frames and the 
thickness of the sheet covering. 


A number- of specimens were tested' in which the geometrical variables 
b and d were systematically varied while R was kept constant and 
equal to l6 inches, The results of these teste yielded a family of ■ 
curves as shown in figure 5^ where the shearing stress in the- sheet cov- 
ering at failure is plotted as a function of the longitudinal spacing b 
for constant values of d/b. An examination of these curves indicates 
that, if the abscissa value of each cxurve . is multiplied by -an appropriate 
expansion factor K, all curves can be made to coincide with, for exam- 
ple, the curve, of d/b equal to O.394. A log-log plot of the expansion 

factor K as a function of d/b indicated that K varies as ^d/b. A 
plot of l/K as a. function of l//d/b - gives the linear variation indi- 
cated in figure 55* A plot of r m»v as a function of is shown in 

figure 56 and indicates that the test results scatter about a common 
curve . 


The next question is> in what manner does the radius R Influence 
the failing load? For the buckling of unstiffened cylinders subjected 
to pure torsion, the experimental results of reference 9 indicate that 
for values of L/R equal to and greater than 3*2 the buckling stress 

Tqj. varies approximately as (t f'R)^ * ^/t^/L as shown in figure 57. It 
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was therefore assumed that^ for identical values of b, d, and Py, 

the critical shearing stress for the general instability of a stiffened 
cylinder would vary as the reciprocal of In order to verify this 

assumption, a number of tests were conducted on l6~inch-radius specimens 
and various values of d and b. A plot of T ynay ^ function of 

/m for the 10- and l6-lnch specimens is shown in figure ^8. Thesi 

results indicated that the assun^jtion was justified and it was therefore 
concluded that the parameter for predicting general inetability in tor^ 
si on is of the form: 


From dimensional reasoning, it follows that the function f(p^,Py) must 

have the dimensions of the J /k power of a length. The simplest assump- 
tion for the function which determines the influence of the section param- 
eters, Pj. and Py, is that it depends on only the geometrical mean 

value y p^ Py. The parameter therefore appears in the form: 


In checking the validity of this parameter, it is necessary to evaluate 
the amount of sheet acting with the frames and longitudinals in order to 
calculate Px and py. It is quite difficult to evaluate by analytic 
means the amount of sheet acting with the reinforcing members; trial cal- 
culations indicated that the best results were obtained if the total 
width of sheet was used. For this reason Px and py were calculated 
with the entire width of sheet assumed to be effective. The variation of 
p^ and Py, with flat sheet assumed, as a function of the effective 

width of sheet Wg is shown in figures 59 to 63 . 

Specimens were also tested in which both the sheet thickness and sec- 
tional properties of the longitudinals were varied. The results of all 
tests are shown plotted in figure 64. It is seen that up to values of 
10,000 povuids per square inch all test values scatter closely about a 
straight line. For higher values of T mpv there is a sudden shift in 
the eap)erimental values. However, the majority of tests again follow a 
straight line having the same slope as the line corresponding to the 
lower values of . Since the observed diagonal-tension field varied 

between about 30 ° to 50 °, it is seen from equation ( 3 ) (see appendix) 
that for a shear stress of 10,000 poxmds per square inch the correspond- 
ing tensile stress would be between 20,000 and 23,000 poimds per square 



f(p^, Py) 
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inch. It was thought that this tensile stress might he sufficiently close 
to the proportional liBiit of the sheet covering to explain the sudden 
shift in the experimental values. However^ an examination of the stress- 
strain curves (figs. 9 to 12) indicates that the tensile stress at a 
shear stress of 10,000 pounds per sq,uare inch is well helow the propoj^- 
tional limit. 

A more desirable presentation of the test data would be to plot 


shearing strain y . Such a presentation would be more general and would 
allow for materials of different physical properties or for changes in 
the physical properties above ihe proportional limit. It had not been re- 
alized at the beginning of the test program that it would be desirable to 
obtain a measure of t/Ge failui’e. For this reason angular deforma- 
tions were measurred on only a rumber of specimens . Not enough measure- 
ments are available to make such a plot. A plot of Tmax/E ^ func— 



in this figure corresponds to that of the sheet covering and was taken as 
lo"^ pounds per s(iuare inch, since this is very close to the actual test 
values obtained for the sheet. 

It should be noted that the curve of figure 65 as presented is 
strictly applicable to the aluminum alloys tested. However, an estimate 
of the failing stress of a reinforced cylinder of different material can 
be obtained by calculating the numerical value of the parameter and as- 
certaining the value of Tjuax/g from the curves of figure 65 . 


The over-all general-instability test program was undertaken to fui^- 
nish the designer with sufficient irJ'ormation to enable him to make an 
estimate of the allowable general-instability stress of a reinforced 
metal cylinder. With the completion of the pure-tors ion- -loading program 
sufficient information is available to estimate the general-instability 
stress for a v-ariety of loading conditions, that is, pure bending# pure 
torsion; combined bending and torsion; and combined bending, traa^ve^-se 
shear, and torsion. 



since Tmax/Gj; corresponds to the 


CONCLUSIONS 
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The quantities involved in the parameter were all varied through a 
sufficient range to verify the generality of the parameter. The sudden 
break which occurs in the summary curve of general- instability failure 
has not been adequately explained. However, it is thought that it is 
connected with a variation in the physical properties of the various 
structural elements of the reinforced cylinder rather than a breakdown 
of the parameter. The reason for this statement is that of the eight 
specimens which failed above a shear stress of 10,000 pounds per square 
inch, one specimen falls on the original curve and the other seven on a 
cvirve which is parallel to the original curve. The parameter as derived 
does not itself give the numerical value of the ultimate stress, but de- 
fines essentially the slope of the ultimate— stress curve. If the shear- 
ing stress at failure t rrinr were replaced by a measured shear-strain 
deformation the experimental results might fall on a single curve. If, 
for example, in the pure— bending parameter *f reference 5 the measured 
unit strain were replaced by a calculated stress corresponding to the 
ultimate bending moment the experimental results would not fall on a 
single curve. 


California Institute of Technology, 

Pasadena, Calif., March 5# 19^6. 

APPENDIX 

INDUCED LONGITUDINAL STRAIN 


In considering a sheet panel bounded by two longitudinal stringers 
and the end mounting rings as shown in figure h(a), it can be seen that, 
upon buckling, the diagonal— tension field produces horizontal components 
of force on the two end rings. The end plate to which one mounting ring 
is bolted is fixed to the base of the testing machine; whereas the other 
end plate is mounted on rollers and is free to move. Therefore, the hor- 
izontal force components must be resisted by the longitudinal stiffeners. 

By considering an element of sheet of unit width (fig.; 4(b)), the 
following relations cein be derived: • 


P = t at cos a 

Py = P sin a = t a^ sin a cos a = t T 


from which 
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= t/sin a cos a = 2T/sin 2a 
= P cos a = t cos^ a = t t cot a 


(3) 


Since the diagonal-tension field exists only beyond the buckling 
stress Tor follows that . 

Pg = t(r - Ter) ® 

The total horizontal component of force around the circumference is 

Pfitot ' ” "^cr^ 

And the unit strain in the longitudinals' is given by: ■■ 

j - 2«t (t - Ter) .cot g- . (^4) 

nA^^ E - , ' ^ • 

where ' • 


n number of longitudinals 
Ag^ area of a longitudinal 


By using the relation of equation (l), equation (4) 
in the form, 


"st 


(Mt - Mlcr) O' 

n-RAgtE 


can be written 
(5) 


In order to compare the values of unit strain given by equation (5) ' 
with experimental values, it is necessary to knpw the value of Mtqj.. 

For purposes of comparison, the value of the moment at which a sharp 
change occurs In the slope of the moment— strain curve, was considered to 
be MTcr* -4 comparison of measured and calculated unit strains Is pre^ 
sented in figures 66 and 67 . For the calculated strains, the value of 
a was assumed to be 45°. It is seen that for the particular stiffeners 
of figure 67 the agreement is quite good; whereas for the two stiffeners 
of figure 66 the agreement is rather poor. ^A check^of available photo- 
graphs indicated that a varied between 30 and 50 depending, on the 
reinforcement spacing. By using the appropriate value of therefore, 

a closer agreement may be obtained. For example, the value of a for 
specimen 208 appears from the photographs to be slightly greater than 
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45 ° and for specimen 209 slightly less than 45°. In both cases, these 
angles tend to improve the agreement between calculated and measured val- 
ues. For specimens I 96 and 200, photographs are not available. 

It should also be noted that there is considerable variation in the 
measured unit strain between longitudinals for a given test speciiaen. 

It is thought that, in general, equation (5) will give values sufficient- 
ly accurate for design purposes provided the correct value of a is used. 
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TABLE I 


Sped— 

men 

Stiff- 

ener 

Frame 

b 

d 

/m 

R 



Px 



/pxPy 

y bd 

- r— — ^3/1 

J PxPy y 

R J 

t 

^max 

max 

1A9 

S2 

> 


r 

5.12 

4 

4.52 

16 

8.0 

36.1 

0.1061 

0.0269 

0.00594 

16.45x10® 

0.010 

10. 70x10" * 

6,600 

162 

S2 




2.56 

4 

3.20 

16 

8.0 

25.6 

.1169 

.0269 

.00648 

25.3 

.010 

13.00 

8,100 

171 

Si 




2.56 

4 

3.20 

16 

8.0 

25.6 

.1145 

.0269 

.00636 

24.82 

.010 

11.85 

7,380 

175 

Si 




2.56 

8 

4.52 

16 

8.0 

36.1 

.1145 

.0234 

.00563 

15.6 

.010 

8.00 

4,970 

176 

Si 




2.56 

2 

2.26 

16 

8.0 

18.075 

.1145 

.0285 

.00670 

37.05 

.010 

15.80 

9,840 

ITT 

Si 




5.12 

16 

8.98 

16 

8.0 

TI.9 

.0967 

.0173£ 

.00373 

5.19 

.015 

5.00 

2,070 

178 

Si 




5.12 

8 

6.36 

16 

8.0 

50.9 

.0967 

.0218 

.00455 

8.94 

.015 

7.50 

3,110 

1T9 

Si 




5.12 

4 

4.52 

16 

8.0 

36.1 

.0967 

.0262 

.00536 

14.8 

.015 

10.90 

4,510 

180 

Si 




5.12 

2 

3.20 

16 

8.0 

25.6 

.0967 

.0289 

.00584 

22.8 

.015 

17.00 

7,040 

181 

Si 




5.12 

16 

8.98 

16 

8.0 

71.9 

.0915 

.01525 

.00318 

4.42 

.020 

7.30 

2,270 

182 

Si 




5.12 

8 

6.36 

16 

8.0 

50.9 

.0915 

.0202 

.00406 

7.97 

.020 

9.25 

2,890 

183 

Si 




5.12 

4 

4.52 

16 

8.0 

36.1 

.0915 

.0251 

.00491 

13.6 

.020 

13.50 

4,190 

184 

Si 




5.12 

2 

3.ao 

16 

8.0 

25.6 

.0915 

.0289 

.00556 

21.7 

.020 

23.35 

7,010 

185 

Si 




5.24 

8 

6.48 

10 

5.62 

37.4 

.0961 

.0218 

.00453 

12.1 

.015 

3.80 

4,025 

186 

Si 




5.24 

4 

4.58 

10 

5.62 

25.8 

.0961 

.0262 

.00533 

20.6 

.015 

6.30 

6,670 

187 

Si 




5.24 

2 

3.24 

10 

5.62 

18.2 

.0961 

.0289 

.00581 

31.9 

.015 

8.90 

9,500 

189 

Si 




10.24 

32 

18.1 

16 

8.0 

144.9 

.0895 

.0143 

.00296 

2.04 

.010 

1.30 

810 

190 

Si 




10.24 

16 

12.8 

16 

8.0 

102.5 

.0899 

.0188 

.00376 

• 3.67 

.010 

2.50 

1,560 

191 

Si 




10.24 

8 

8.98 

16 

8.0 

71.9 

.0899 

.0234 

.00455 

6.33 

.010 

4.10 

2,550 

192 

Si 




10.24 

4 

6.36 

16 

8.0 

50.9 

.0899 

.0269 

.00514 

10.1 

.010 

6.76 

4,210 

193 

Si 




2.56 

2 

2.26 

16 

8.0 

18.075 

.1145 

.0285 

.0067 

37.05 

.010 

16.00 

9,940 

194 

S3 




5.24 

8 

6.48 

10 

5.62 

37.4 

.0947 

.0228 

.00466 

12.45 

.0115 

2.90 

4,010 
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TABLE I - Coaocludad- 


Speci- 

men 

Stiff- 

ener 

Frame 

b 

d 


E 



Px 

Py 


/p P 

y M 

r /p P~ 1 
R 


t 

*^max 


195 

s. 




5.24 

4 

4.58 

10 

5.62 

25.8 

0.0947 

0.0266 

0.00533 

20.65x10® 

0.0115 

4 . 95 X 15 ’* 

6,780 

196 

S 3 




5.84 

2 

3.24 

10 

5.62 

18.2 

.0947 

.0287 

.00571 

31.30 

.0115 

6.70 

9,270 

197 

Sz 




2.62 

8 

4.58 

10 

5.62 

25.8 

.1166 

.0228 

.00558 

21.60 

k 

.0115 

5.10 

7,050 

198 

Sa 




2.62 

4 

3.24 

10 

5.62 

18.2 
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Figure 13.- Specimen with deflection pointers. 


1 



NACA TN No. 1197 Fig. 14 



Failure of a specimen with 2-inch 
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FIGURE 19. 
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Specimen No. 198 

FIGURE 47 
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SHEAR STRESS “T X I0‘^ 
FIGURE 49. 
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FIGURE 64- GENERAL INSTABILITY FAILURE; PURE TORSION. 
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